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SUMMARY 



Two components of a structure which are located side by side, 
will come in contact by certain force and will transfer the 
compressive force along the contact area. If the force acts in the 
opposite direction, the elements will separate and no force will be 
transferred. If this contact is modelled, the load path will be 
correctly represented, and the load redistribution results in more 
realistic stresses in the structure. This is accomplished by using 
different sets of rigid elements for different loading conditions, or 
by creating multipoint constraint sets. Comparison of these two 
procedures is presented for a 4 panel unit (PU) stowage drawer 
installed in an experiment rack in the Spacelab Life Sciences (SLS- 
2) payload. 


INTRODUCTION 

The Spacelab is a reusable laboratory that is carried in the 
cargo bay of the Space Shuttle. Experiments in several different 
disciplines such as astronomy, life sciences, and material science 
are accommodated in this modular laboratory for various Shuttle 
missions. The experiment hardware is mounted in the experiment 
racks located in either side of the module, in overhead lockers, or in 
the center aisle, as shown in Figure 1. 

4PU STOWAGE DRAWER 
Configuration 

The 4 Panel Unit Stowage Drawer is mounted in the experiment 
rack used in the SLS-2 Mission. The experiment equipment and the 
accessories are stowed in the drawer. The finite element model of 
the drawer, with its coordinate system, is shown in Figure 2. The 
drawer is connected to the slide with 6 screws on each slide, and 
the slides are connected to the rack posts. The front panel is latched 
to the front rack posts. Two configurations of the slides are 
examined. 
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Case 1: The contact surface is normal to the X-axis, as shown 
in Figure 4, which is the actual configuration. The slide shown is 
schematic, and not the actual slide. 

Case 2: The contact surface is inclined. This is achieved by 
raising the slide by 12.7 mm as shown in Figure 5. 

The Method of Modelling the Contact 

During liftoff and landing flight events, the Shuttle and its payload 
are exposed to quasi-static and random loads. The +X force brings 
the right slide and drawer in contact. As a result, this force is 
transferred to the slide throughout the length of the slide and not 
just by the screws. When the force acts to the left (-X), the contact 
along the length is lost and the right slide is connected by screws 
only. This time the contact takes place between the left slide and 
the drawer. 


Generally, this is modelled using rigid elements. For load case 
101, which includes +X force (see Table 1), ail the contacts are 
modelled between the right slide and the drawer, and the analysis is 
completed. For the load case 103, which includes -X force, contact 
will be modelled between the left slide and the drawer with a new 
set of rigid elements, removing the old set of elements, and a second 
analysis will be performed. This means post-processing will be 
performed on two output files. The rigid elements simulating 
contact are shown in Figure 3. These are included in the analyses, as 
needed. 

Alternately, the contact is modelled with a multipoint 
constraint equation in place of the rigid element. In this method, a 
different set of MPC equation can be written for a different subcase, 
resulting in a single analysis for multiple subcases. 
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CASE 1 . CONTACT SURFACE NORMAL TO THE GLOBAL X-AXIS 


Modelling of Contact by MPC Equation 

A rigid link is used to write the MPC equation, as shown in 
Figure 6 (ref 1.). Since the physical contact cannot resist moments, 
no rotations will be allowed at the end of the rigid link. In the 
current case, the link is horizontal, i.e., AL - AX ■ u-j . 

The MPC equation is uia* u 1B ■ 0. The MPC set 1 is written for the 

subcase 101 to represent right slide contact, and the MPC set 2 is 
written for the subcase 103 to represent the left slide contact. The 
MPC equations and the MPC forces are shown in Table 2. The grid 
points shown are on the slide. Grid points 471 and 472 show forces 
in opposite directions, indicating tension and lack of contact. In this 
situation, these equations should be removed and reanalysis must be 
performed. In the current analysis, this is not pursued. 

Modelling with CRIGD2 Elements 
The elements modelled and the results for subcase 101 and 
subcase 103 are shown on Tables 3 and 4, respectively. The 
dependent degree of freedom is 1. The equation generated 
corresponds to row 1 of equation 56 (ref. 2) shown below. 
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( 56 ) 


In the equation, (Zb - Za) corresponds to Zab on Figure 6 which is 
zero for this case, and (Yb -Ya) is also zero. Hence, the equations 
generated are the same as the MPC equations and the results from 
both the analyses will be identical. 

Discussions of the Two Methods 
As expected, the results from both of these methods are the 

same. 
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CASE 2. CONTACT SURFACE INCLINED TO THE GLOBAL X-AXIS 


Only subcase 101, which involves +X and +Z loads, will be used 
in the following analyses. Due to these forces, contact will be made 
in the X and Z directions as shown in Figure 5. 

Modelling wit h MPC Equations 

MPC equations are written to satisfy the geometry of the rigid 
links shown in Figure 6 . As stated before, no rotation will be 
allowed. 

Xab z AB 

Hence AL - -p- ui + U 3 (1) 

from the geometry of the inclined link in Figure 5, 

Xab - 18.606 mm, Zab - 12 7 mm . L * 22 -527 mm 
Substituting in equation (1) 

AL « .8259 u-| + .5638 U3 
ALa ■ .8259 u-)A + -5638 U 3 A 
ALp- .8259 u-jb + -5638 U 3 B 

to satisfy the condition ALa - ALb - 0, the MPC equation is 
.8259 (uia -uiB) + -5638 (U 3 A - U 3 B ) * 0 ( 2 ) 

Table 5 shows all the MPC equations input for all the contacts, 
followed by the forces of multipoint constraint, in the grid points on 
the slides. 


Modelling w ith CRIGD2 Elements, 

CRIGD2 are modelled with components 1 and 3 as dependent 
degrees of freedom to simulate the contact in X and Z directions. The 
constraint equations generated correspond to rows 1 and 3, in the 
equation 56. The term (Zb -Za) in row 1 and (Xb -Xa) in row 3 are 
non-zero. These terms correspond to component 5, and it is expected 
that constraint moment R2 will be generated. The list of elements 
and the results are tabulated in Table 6. 
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Modelling w ith CRIGDR Elements 


CRIGDR elements are modelled with component 1 as the 
dependent degree of freedom. The remaining 5 translational 
components are considered as reference degrees of freedom (ref. 2). 
Equation 48 (ref. 2) is used in the element formulation shown below. 

(UA1 - U B1 ) *1 + ( U A2 * U B2) ! 2+ ( U A3 -U B3 ) ! 3 " 0 ( 48 ) 

XAB , ZAB 

In this equation, direction cosine I2 - 0, li - , I3 - 

which essentially is MPC equation (2), and the results from this 
analysis will be same as from the MPC equation. 

A list of the elements and the results are tabulated in Table 7. 

Comparison of the T hree Analyses 

It is shown that the formulation of MPC equations and the 
CRIGDR are identical, and the results tabulated in Tables 5 and 7 are 
identical as expected. The CRIGD2 results are different than the 
other two because this involves rotations. In this instance, R2 
moments are generated as expected and the Z components are off by 
about ±20 percent. 


CONCLUSIONS 

The best way to model contact is by writing MPC equations 
since a single analysis, is possible for multiple subcases. CRIGDR is 
the second choice. 
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TABLE i TOTAL APPLIED POBCE OH THE STHUCTUBE ( HEWTOHS I 
DZIECTZOH XT Z 

SUBCASE 101 1001.2 175.3 27S.4 

SUBCASE 103 -117.0 175.3 27S.4 

TABLE 2 PABTIAL IHPUT AHD BCSULTS - CASE 1 
SUBCASES 101.103 - USE OP KPC EQUATZOHS 
$ KPC EQUATZOHS TO SIMULATE COHTACT ZH X-DZB 
$ POB SUBCASE 101 CASE 1 


NFC 


1 


454 1 

1.0 

1070 


1 

-1 

.0 

MFC 


l 


471 1 

1.0 

1104 


1 

-1 

.0 

MFC 


1 


472 1 

1.0 

1 13*8 


1 

-1 

.0 

MFC 


1 


473 1 

1.0 

1172 


1 

-1 

.0 

MFC 


1 


474 1 

1.0 

1204 


1 

-1 

.0 

MFC 


1 


475 1 

1.0 

1240 


1 

-1 

.0 

MFC 


1 


474 1 

1.0 

1251 


1 

-1 

.0 

MFC 


1 


477 1 

1.0 

1342 


1 

-l 

.0 

MFC 


1 


471 l 

1.0 

1375 


1 

-l 

.0 

$ MFC EQUATIONS TO 

i SIMULATE COHTACT ZH -X 

-DIR 





$ POt SUBCASE 

103 

CASE a 







MFC 


2 


5454 1 

1.0 

5070 


1 

-1 

.0 

MFC 


2 


5471 1 

1.0 

6104 


1 

-1 

.0 

MFC 


2 


5472 1 

1.0 

6138 


1 

-1 

.0 

MFC 


2 


5473 1 

1.0 

6172 


1 

-1 

.0 

MFC 


2 


5474 1 

1.0 

6206 


1 

-1 

.0 

MFC 


2 


5475 1 

1.0 

6240 


1 

-1 

.0 

MFC 


2 


5474 1 

1.0 

6251 


1 

- 1 . 

ft 

MFC 


2 


5477 1 

1.0 

6342 


1 

-1 

.0 

MFC 


2 


5471 1 

1.0 

6376 


1 

-1 

.0 

F O 

BCE 

S 

o r 

M U L T I - F 

O X N T 

C O 

N S 

T ! 

5 A I 

N T 

*t?*C**E 

101 


CASK 1 







POINT ID 

l. TTFE 

T1 

T2 

T3 

51 

52 

53 


471 


a 

2 . 000543E+00 

0.0 

0.0 

0.1 


0.0 

0.0 


472 


G 

5 . *274731-01 

0.0 

0.0 

0.4 


0.0 

0.0 


473 


G 

-2 . 36 I567E-0 1 

0.0 

0.0 

0.4 


0.0 

0.0 


474 


G 

-4.407188E-01 

0.0 

0.0 

0.4 


0.0 

0.0 


475 


G 

-4 . 436086E+00 

0.0 

0.0 

0.4 


0.0 

0.0 


474 


G 

-5.1991 1 2E+00 

0.0 

0.0 

0.4 


0.0 

0.0 


477 


G 

-2 .6444221*00 

0.0 

0.0 

0.4 


0.0 

0.0 


4 71 


G 

-4 . 0577601*00 

0.0 

0.0 

0.4 


0.0 

0.0 

SUBCASE 

103 


CASE 1 







FOIST ZO. 

TYPE 

T1 

T2 

T3 

51 


52 

53 


5471 


G 

-1 .5055141+00 

0.0 

0.0 

O.C 


0.0 

0.0 


5472 


G 

-7.575473E-01 

0.0 

0.0 

o.c 


0.0 

0.0 


5473 


G 

3 . 0011 1 2E-0 1 

0.0 

0.0 

o.c 


0.0 

0.0 


5474 


G 

1 .1035401-01 

0.0 

0.0 

o.c 


0.0 

0.0 


5475 


G 

1 .7912771.0® 

0.0 

0.0 

o.c 


0.0 

0.0 


5475 


G 

4. B3SSS7E.00 

0.0 

0.0 

o.c 


0.0 

0.0 


5477 


G 

2 . 1519O9E.O0 

0.0 

0.0 

o.c 


0.0 

0.0 


5471 


G 

3.2I3999E.00 

0.0 

0.0 

o.c 


0.0 

0.0 
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TABLE 3 PAXTXAL IBPUT ABO BESULTS - CASE 1 

SUBCASE 101 - USE OP CXXCD2 ELEMESTS • 

$ BXSXD ELENEBTS MODELED TO SIMULATE COBTACT XB X-DIX 
$ POX LOADCASE 101 CASE 1 


CRXG02 

410 

454 

1070 

1 





C1XOD2 

411 

471 

1104 

1 





CXXGD2 

442 

472 

1134 

1 





C&XGD2 

443 

473 

1172 

1 





C1ZGD2 

444 

474 

1206 

1 





C1XG02 

445 

475 

1240 

1 





CRIGD2 

446 

476 

1291 

1 





C1XOD2 

447 

477 

1342 

1 





CRIGD2 

444 

474 

1376 

1 





SUBCASE 101 



CAS El 






P 0 R C t 5 

O P 

M U L 

T X - P 

0 I H T 

c 

OHS 

T t A 

I W T 

POIITP ZD. 

TYPE 


T 1 

T2 

T3 

R1 

R2 

R3 

471 

G 

2 . 0006 6 3C+00 

0.0 

0.0 

0.0 

0 . 0 

0.0 

472 

G 

9. 6 27473E -01 

0.0 

0.0 

0.0 

0.0 

0.0 

473 

G 

-7O69567E-01 

0.0 

0.0 

0.0 

0.0 

0.0 

474 

G 

-4.4071 8 IE-0 1 

0.0 

0.0 

0.0 

0.0 

0.0 

473 

G 

-4 . 4 3 60S 6E+00 

0.0 

9.0 

0.0 

0.0 

0.0 

476 

G 

-5 . 499 11 2E+00 

0.0 

9.0 

0.0 

0.0 

0.0 

477 

G 

-2.644422E+00 

0.0 

9.0 

0.0 

0.0 

0.0 

471 

G 

-4 . 057767000 

0.0 

9.0 

0.0 

0.0 

0.0 

TA1LE 4 

PARTIAL 

IBPUT 

AH 0 RESULTS -CASE 1 




SUBCASE 103 - 

USE OP 

CRXGD2 

ELCMEHTS 






S XXGXD ELEMEBTS MODELED TO SIMULATE COBTACT IB -X-DIX 


S POX SUBCASE CASE 103 CASE 1 


CXXGD2 

5440 

5454 

6070 


1 





CXXGD2 

5441 

5471 

6104 


1 





C1XGD2 

5442 

5472 

6134 


1 





CXXGD2 

5443 

5473 

6172 


1 





CXXGD2 

5444 

5474 

6206 


1 





CRXGD2 

5445 

5475 

6240 


1 





CXXGD2 

5446 

5476 

6291 


1 





CHXGD2 

5447 

5477 

6342 


1 





CBXGD2 

5444 

5474 

6376 


1 





P O R C E 

S OP 

H U L T I - P 

O * H 

T C 

O If s 

T R A 

I w 

POXHT ID. 

TYPE 

T1 

T2 


T3 

R1 

R2 

R3 

5471 

G 

-I . «0»S14E.OO 

0.0 


0.0 

0.0 

0 . 0 

0.0 

5472 

G 

-7.573473E-01 

0.0 


0.0 

0.0 

0.0 

0.0 

5473 

G 

3 .0041121-01 

0.0 


0.0 

0.0 

0 . 0 

0.0 

5474 

G 

1.1033*01-01 

0.0 


0.0 

0.0 

0.0 

0.0 

5475 

G 

3.7M277E+00 

0.0 


0.0 

0.0 

0.0 

0.0 

5476 

G 

4 . S 3 1* *71+00 

0.0 


0.0 

0.0 

0.0 

0.0 

5477 

G 

2 .151* 09E+00 

0.0 


0.0 

0.0 

0.0 

0.0 

5474 

G 

3 . 2I3PP9E.00 

0.0 


0.0 

0.0 

0.0 

0.0 
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$ MFC 

$ FOR 

MFC 

♦MFC1 

MFC 

♦MFC2 

MFC 

♦MFC3 

MFC 

♦MFC 4 

MFC 

♦MFCS 

MFC 

♦MFCS 

MFC 

♦MFC? 

MFC 

♦MFCS 

MFC 

♦MFCS 


TABLE 5 PARTIAL INPUT AMD RESULTS -CASE 2 




SUBCASE 

ioi - use or MFC 

EQUATIONS 


EQUATIONS TO SIMULATE CONTACT IN X 

AND Z DIR 


SUBCASE 

101 

CASE 2 

(INCLINED SURFACE) 



1 

454 

1 0.1259 

1070 1 -0.4259 

♦ MFC1 



454 

3 0.5434 

1070 3 -0.5434 



1 

471 

1 0.4259 

1104 1 -0.4259 

♦MFC2 



471 

3 0.5434 

1104 3 -0.5634 



1 

472 

1 0.4259 

1134 1 -0.4259 

♦MFC3 



472 

3 0.5434 

1134 3 -0.5634 



1 

473 

1 0.4259 

1172 1 -0.4259 

♦MFC4 



473 

3 0.5434 

1172 3 -0.5634 



1 

474 

1 0.4259 

1206 1 -0.4259 

♦MFC5 



474 

3 0.5434 

1206 3 -0.5634 



1 

475 

1 0.4259 

1240 1 -0.1259 

♦MFCS 



475 

3 0.5434 

1240 3 -0.5634 



1 

474 

1 0.4259 

1291 1 -0.4259 

♦MFC7 



474 

3 0.5434 

1291 3 -0.5634 



1 

477 

1 0.4259 

1342 1 -0.4259 

♦MFCS 



477 

3 0.5434 

1342 3 -0.5434 



1 

471 

1 0.1259 

1376 1 -0.4259 

♦MFCS 



471 

3 0.5434 

1376 3 -0.5634 


SUBCASE 

101 

CASE 2 (INCLINED SURFACE ) 


r o 

R C 

E S 0 

F MULTI- 

FOIST COKSTftA 

I N T 

POINT ID 


TTFE 

T1 

T2 TJ K1 

R2 ] 

471 


G 

2 . 34 94 4 2E + 00 

0.0 1.CO4145E.O0 0.0 

0.0 0 

472 


G 

1.092 4 54E^0 0 

0.0 7 . 4t 0 3 S7C-0 1 0.0 

0.0 0 

473 


0 

1 . 3 0 730 3E + 00 

0.0 t . 52 4 29 7E-0 1 0.0 

0.0 0 

474 


G 

-1.05 1737E+00 

0.0 -7.179473E-01 0.0 

0.0 0 

475 


G 

—4 . 54292SE+00 

0.0 -3 . 1012241+00 0.0 

0.0 0 

4 74 


G 

-4 . 415951E^OO 

0.0 -4 .379I44E.00 0.0 

0.0 0 

477 


G 

-2 . 20 24 2 7E + 00 

0.0 -1 . 503 4I5E.0 0 0.0 

0.0 0 

471 


G 

-3 . 1 4 7445E^OO 

0.0 -2.34tC15E.00 0.0 

0.0 0 , 
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TABLE 6 PARTIAL INPUT AND RESULTS - CASE 2 
SUBCASE XOl - USE OF CRXGD2 ELEMENTS 


$ RIGID ELEMENTS 
$ DUE SUBCASE 101 

MODELED TO 
CASE 2 

CONTACT 

IS Z AMD Z- 

CRXGD2 

410 

454 

1070 

13 

CRZGD2 

4 ft 1 

471 

1104 

13 

CRIGD2 

412 

472 

1138 

13 

CRIGD2 

413 

473 

1172 

13 

CRXGD2 

4 ft 4 

474 

1206 

13 

CRXGD2 

415 

475 

1240 

13 

CRXGD2 

486 

476 

1291 

13 

CRIGD2 

487 

477 

1342 

13 

CRXGD2 

488 

478 

1376 

13 

SUBCASE 

101 

CASE 

2 ( I8CLIBED SURFACE ) 


FORCES 

O F 

M U L T I - P 

O I B T 

C O B S T R 

A I H 

T 


POINT XD. 

TTPE 

T1 

T2 

T3 

R1 

12 

R3 

471 

G 

2 . 494523E+00 

0.0 

2 • 1218 2ftf +00 

0.0 

3 . 06 9774C— 01 

0.0 

472 

G 

1.101543E+00 

0.0 

9 . 2 8164 4E-01 

0.0 

1 . 29 101 7E-0 1 

0.0 

473 

G 

1.3203 3 7E+00 

0.0 

5.0131 ME-fci 

u.o 

-2.929S34C-01 

0.0 

474 

G 

-1.060023E+00 

0.0 

- 5 . 1 1 4 4 7 IE-0 1 

0.0 

1 .040M2E-01 

0.0 

475 

G 

-4.53445 5 E + 00 

0.0 

-2 . 13191 1 E+00 

0.0 

1.I77255E“01 

0.0 

476 

G 

-6 .43764 3E+00 

0.0 

-3.74 24131+00 

0.0 

4 . 77503 IE-0 1 

0.0 

477 

G 

- 2 • 2 205 1 2E + 00 

0.0 

-1 . 2ft 2993E + 00 

0.0 

1.704 6 3 5E-0 1 

0.0 

478 

G 

“3 .182428 E + 00 

0.0 

“1 • 6709 8 1E+00 

0.0 

3 . 6722 0 5E-0 1 

0.0 

TABLE 7 

PARTIAL 

INPUT AND RESULTS CASE 2 




SUBCASE 101 -USE O T CR1CB* ELEKEBTS 





$ RIGID ELEMERTS TO 

SIMULATE CONTACT IN X 

AND Z -DIR 




$ FOR SUBCASE 101 

CASE 2 (INCLINED SURFACE) 




CRXGDR 

480 

454 1070 

1 





CRXGD1 

481 

<71 1104 

1 





CRXGDR 

482 

472 1131 

1 





CRXGDR 

483 

473 1172 

1 





CRXGDR 

484 

474 1 2 0< 

1 





CRXGDR 

415 

475 1240 

1 





CRXGDR 

486 

47< 1291 

1 





CRXGDR 

487 

477 1342 

1 





CRXGDR 

488 

471 137C 

X 





SUBCASE 101 


CASE 2 ( IBCLIBED P LABE ) 




F O R C E S 

O F 

M U L T I - P 

O I B T 

C O B S T R 

A I 8 

T 


POXBT ID. 

TTPE 

T1 

T2 

T3 

XI 

R2 R3 


471 

G 

2 . 37IK3E+00 

0.0 

1 . 6 23320E + 00 

0.0 

0.0 0.0 


472 

G 

1 . 119 243E + 00 

0.0 

7 . (39 II 3E-01 

0.0 

0.0 0.0 


473 

G 

1 .316421E+00 

0.0 

1 .9I5I10E-01 

0.0 

0.0 0.0 


474 

G 

-1 .05334 7E+00 

0.0 

-7. 190014 E-01 

0.0 

0.0 0.0 


475 

G 

—4 . 54 3 00 (E+00 

0.0 

-3 . 10 10 2 1 E + 00 

0.0 

0.0 0.0 


476 

G 

“6 . 414644E+00 

0.0 

-4 . 3 7ft 59 6 E + 00 

0.0 

0.0 0.0 


477 

G 

-2 . 1 83 3 3 7E + 00 

0.0 

-1 . 490 3 3 3E + 00 

0.0 

0.0 0.0 


478 

G 

“3 . 126624E + 00 

0.0 

-2 . 13 4214E + 00 

0.0 

0.0 0.0 
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FU STOWAGE DRAW A ASSEMBLY 
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FIG. 2 FINITE ELEMENT MODEL 




FIG 3. PLAN VIEW SHOWING RIGID ELEMENTS 
CONNECTING DRAWER TO SLIDE 
( The numbers shown are grid points on the slide) 
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